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Abstract 

We study a model for the pinning of vortices in a two-dimensional, inhomogeneous, Type-II 
superconductor in its mixed state. The model is based on a Ginzburg-Landau (GL) free energy 
functional whose coefficients are determined by the mean field transition temperature T c q and 
the zero temperature penetration depth A(0). We find that if (i) T c q and A(0) are functions of 
position, and (ii) A 2 (0) oc T c y with y > 0, then vortices tend to be pinned by regions where T c o, 
and therefore the magnitude of the superconducting order parameter A, are large. This behavior 
is in contrast to the usual picture of pinning in Type-II superconductors, where pinning occurs in 
the small-gap regions. We also compute the local density of states of a model BCS Hamiltonian 
with d-wave symmetry, in which the pairing field A is obtained from the Monte Carlo simulations 
of a GL free energy. Several features observed in scanning tunneling spectroscopy measurements 
on YBa2Cu30g+x and Bi2Sr2CaCu208+:r are well reproduced by our model: far from vortex cores, 
the local density of states (LDOS) spectra has a small gap and sharp coherence peaks, while near 
the vortex cores it has a larger gap with low, broad peaks. Additionally, also in agreement with 
experiment, the spectra near the core does not exhibit a zero energy peak which is, however, 
observed in other theoretical studies. 
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I. INTRODUCTION 



It is generally believed that vortices in Type II superconductors tend to be pinned in 
regions where the gap is small [1!]. This is true because a vortex, being a region where the 
gradient of the superconducting order parameter is large, locally increases the gradient part 
of the free energy. Since this local increase is less in regions where the gap is smaller than 
average, the vortex tends to migrate to such regions, according to this picture. 

In this paper, we describe a simple model for vortex pinning in an inhomogeneous two- 
dimensional (2D) superconductor, in which the vortices tend to be pinned in regions where 
the gap is larger than its spatial average. The model is based on a Ginzburg-Landau (GL) 
free energy functional in which both the mean-field transition temperature T c0 and the 
zero-temperature penetration depth A(0) are functions of position, but are correlated in 
such a way that regions with large T c0 also have large A(0). This assumption seems to 
apply to some of the high-T c cuprate superconductors: according to scanning tunneling 
microscopy (STM) experiments on cuprates 
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, |9J], regions that have a 

large gap (proportional to T c0 in this model) also have a single-particle local density of 
states (LDOS) with low, broad peaks, suggestive of a low superfluid density in these regions 
[proportional to 1/A 2 (0)]. Our model is a generalization of an earlier approach intended to 
treat inhomogeneous superconductors in zero magnetic field [lo|. 

In order to further test this vortex pinning model, we also examine the quasiparticle 
LDOS near the vortex cores within this model. The proper theoretical description of this 
LDOS near the cores is one of the unsolved issues in the field of high-T c superconductivity. 
STM experiments on YBa 2 Cu 3 6+x (YBCO) fl [l2| and Bi 2 Sr 2 CaCu 2 8+;c (Bi2212) 13, 
14] . 15 . 16], show that the local density of states (LDOS) near vortex cores in those materials 



has the following characteristics: a dip at zero energy, small peaks at energies smaller than 
the superconducting gap, and low, broad peaks at energies above the superconducting gap. 
Those features contrast with the spectra shown by conventional superconductors near vortex 
cores, where the LDOS usually has a peak at zero energy in clean superconductors (i.e., 
those with a mean- free path larger than the coherence length) [l?], [3] , or is nearly energy- 
independent in dirty superconductors [la. Il9|. 

The zero-energy peak in the LDOS near vortex cores of clean conventional superconduc- 
tors can be understood in terms of electronic states with subgap energies bound to vortex 
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21] . However, the structure of the spectra near vortex cores of cuprates is still 



lacking an explanation. Several authors have suggested that this structure is due to some 
type of competing order which emerges within the vortex cores when superconductivity is 



suppressed by a magnetic field 
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script ions of the spectra near vortex cores, have used t 



to solve various microscopic Hamiltonians 
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261 ] . Some of those models, and other de- 
le Bogoliubov-de Gennes method 
3JJ , such as BCS-like models 
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with d-wave symmetry. One recent model for the LDOS near the vortex cores, proposed 
by Melikyan and Tesanovic [3o| , used a Bogoliubov-de Gennes approach to a tight-binding 
Hamiltonian. These authors found that, if a homogeneous pairing field is used in the mi- 
croscopic Hamiltonian, the LDOS exhibits a zero-energy peak on the atomic sites that are 
closest to the vortex cores; this peak is, however, absent from all other sites near the vortex 
cores. On the other hand, they found that introducing an enhanced pairing strength for 
electrons on nearest-neighbor atomic sites near the vortex cores leads to a suppression of 
the zero energy peak, thus obtaining a better agreement with experiment. They speculate 
that the enhancement could be due to impurity atoms that pin the vortices, to a distortion 
of the atomic lattice by the vortex itself, or to quantum fluctuations of the superconducting 
order parameter. The assumption of an enhanced pairing strength is consistent with the 
model that we describe in this paper. 

In our model, having introduced this correlation between the GL parameters, we anneal 
the system to find both the magnitude and the phase of the superconducting order pa- 
rameter that minimize the GL free energy at low temperatures. Since the magnetic vector 
potential enters the GL free energy functional, this procedure naturally leads to vortex for- 
mation. The vortex cores can be identified in our simulations as regions with a large phase 
gradient. We find that, in inhomogeneous systems, vortices tend to be pinned in regions 
where the superconducting gap is large. For comparison, we also perform a similar annealing 
procedure for homogeneous systems. In this case, contrary to inhomogeneous systems, the 
magnitude of the superconducting order parameter is reduced near vortex cores. Thus, the 
assumption that the gap is large in regions with small superfluid density in inhomogeneous 



systems, originally intended to model superconductors in a zero magnetic field [10j, leads 
naturally to pinning of the vortices in large-gap regions. This approach might therefore be 
complementary to that of Ref. [30] mentioned above. 

To connect our vortex pinning model to previous studies of the LDOS near vortex cores, 
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we have also studied a microscopic Hamiltonian for electrons on a lattice. This is a tight- 
binding model in which electrons on nearest neighbor sites experience a pairing interaction 
of the BCS type with d-wave symmetry l(J,|32j|. The LDOS is obtained by exact numerical 
diagonalization of this Hamiltonian. We take the pairing stregth between electrons on near- 
est neighbor sites to be proportional to the value of the superconducting order parameter 
as determined by the GL simulations using the GL functional just described. 

Using this combination of a GL free energy functional and a microscopic d-wave BCS 
Hamiltonian, we find that a number of experimental results on cuprates are reproduced well 
by our model for inhomogeneous systems. For example, far from vortex cores, the LDOS 
shows sharp coherence peaks. Near the vortex cores, our calculated LDOS does not show 
a spurious zero energy peak; instead it exhibits a large gap, as well as low, broad peaks 
which occur at energies larger than the value of the superconducting energy gap observed 
far from vortex cores. Also, in agreement with experiment, the LDOS curves near vortex 
cores are similar to those in the large gap regions in systems with quenched disorder but zero 
magnetic field llOj); this connection is discussed in section UTTl One feature not captured by 
our model is the existence of small, low energy peaks in the LDOS observed near the vortex 
cores. 

The rest of the present article is organized as follows. In Section [Til we present the 
GL free energy functional, as well the microscopic Hamiltonian. In Section IHIt we present 
results for homogeneous and inhomogeneous systems in a magnetic field, as well as results 
for inhomogeneous systems in a zero magnetic field for comparison. Finally, in Section IIVI 
we conclude with a discussion and a summary of our work. 
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II. MODEL 



A. Ginzburg-Landau free energy functional 

We use a model for a single layer of a cuprate superconductor in a perpendicular magnetic 



field based on a GL free energy functional of the form described previously! 101]: 
ST = 5 (u + 3 ) W2 + 2(9.38) A?(0)t^ W " 

In Eq. ([1]) the first and second sums are carried over M square cells, each of area equal to 
the zero-temperature GL coherence length squared, £q, °f a square lattice into which the 
superconductor has been discretized for computational purposes. The third sum is carried 
out over nearest neighbor cells (ij). Here 

1>i - (2) 

fro 

where Aj is the complex superconducting order parameter of the i th cell, Eq is an arbitrary 
energy scale (which we take to be the hopping constant th op ~ 200 meV) 

t=*-f, (3) 
fro 

is the reduced temperature, T is the temperature, and ks is the Boltzmann constant. Also, 
Aj(0) is the T = penetration depth and t c0i = kBT c0i /E , where T c0i is the mean-field 
transition temperature of the i th cell. Therefore, in discretizing the superconductor, we have 
assumed that A(0), T c0 and A are constant over distances of order £ . Finally, 

% = f c \l A ^ • d? ^ 

is the integral of the vector potential A(r) from the center of cell i, located at i, to the center 
of cell j, located at j, and K x = ^ 4 rf/[32(9.38)7rm*V|], where fi% ~ 5.4 x 10~ 5 eV-A 3 is the 
square of the Bohr magneton, m* is twice the mass of a free electron, e is the absolute value 
of its charge, and d is the thickness of the superconducting layer. In all of our simulations 
we use periodic boundary conditions. A gauge for the magnetic field that allows this is given 



in Ref. 

N v = l] 
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33j [see, e. g., eq. (51) of Ref. [10|, which gives this gauge choice explicitly for 



We note that the connection between the local superfluid density n s ^(T) and penetration 



depth Aj(T) is such that n g) i(0) oc 1 / Af (0) [10|. Thus, the quantity 1/A 2 in eq. (1) is a way 



of describing the local superfluid density, which varies over a length scale of £ - 

We will be using the above free energy functional at both T = and finite T. Although 
we call this a "Ginzburg-Landau free energy functional," this name is really a misnomer, 
since the original GL functional was intended to be applicable only near the mean-field 
transition temperature. Strictly speaking, the correct free energy functional near T = 
should not have the GL form, but would be expected to contain additional terms, such as 
higher powers of \ip\ 2 - We use the GL form for convenience, and because we expect it will 
exhibit the qualitative behavior, such as vortex pinning in the large-gap regions, that would 
be seen in a more accurate functional. 

We assume that the magnetic field is uniform, which is a good approximation for cuprate 
superconductors in their mixed state, provided the external field is not too close to the lower 
critical field H c \. In the cuprates, the approximation is satisfactory because the penetration 
depth is of the order of thousands of A, while the intervortex distance for the fields we 
consider is ~ 100 A. We employ a gauge that permits periodic boundary conditions, such 



that the flux through the lattice can take any integer multiple of hc/e lOj, l33| . Thus, the 
number N v of flux quanta hc/(2e) must be an integer multiple of two. 



The procedure for choosing the parameters Aj(0) and t c oi is similar to that used in Ref. 10] . 
Basically, in most of the system (a regions), we take Aj(0) ~ A(0), where A(0) is the in- 
plane penetration depth of a bulk cuprate superconductor [A(0) ~ 1800 A in LM2212, for 
example], while t c oi is determined from the typical energy gap in the LDOS as observed in 
STM experiments. However, we will also introduce regions {(3 regions) in which t cQ and A(0) 
are larger than those bulk values. Throughout this article, when we refer to a "gap in the 
LDOS," we mean the distance between the two peaks in the LDOS spectra. This is the 
same gap definition used in Ref. {4]. 

In Ref. 10j, we generally introduced inhomogeneities in the superconducting order pa- 
rameter Aj by assuming a binary distribution of t ci0 , randomly distributed in space: a cells 
with a small t ci0 and f5 cells with a large i ci0 . We also assumed a correlation between Aj(0) 



and t c0i of the form 

A 2 (0) = A 2 (0)^Y, (5) 



with y = 1. More generally, Eq. ([5]), with y > 0, accounts for the fact that in STM 



experiments, regions with a large gap seem to have a small superfluid density (low and 
broad peaks.) 

In the present article, instead of randomly distributed (3 cells, we introduce two square 
regions with only (3 cells, while the rest of the lattice is assumed to have only a cells. We 
make this choice to study the effects of these inhomogeneities on field-induced vortices. We 
also consider a more general model than Ref . [h| , allowing y to have values other than only 

y = i. 

As shown in Eq. (36) of Ref. [10|, in the absence of thermal fluctuations the coupling 
constant JxY,ij between cells i and j is approximately 

T M 2(9.38) v /(l-t/t e0a )(l-t/t c0j ) ^ 

JxY «® * xMW) Kl - (6) 



(The factor of K% is missing in Ref. }10|.) If t « t c0i and t « t c0 j, then 

, , 2(9.38)^1 

If we further choose a binary distribution of t c0i , that is, 

{t c o, if i is on an a cell, 
(8) 
ftco, if i is on a (3 cell, 

where / is any positive number (typically / > 1), then 

1 if i and j G a, 



2(9.38)^ 



-V if « € a and j G /3 or if i G /3 and j G a, (9) 
if z and j G /3. 

This expression shows that in regions with a large gap the coupling between XY cells is 
small, reflecting the large penetration depth in those regions. 

B. Microscopic Hamiltonian 

Besides using a GL free energy to explore vortex pinning in an inhomogeneous supercon- 
ductor, we have also studied the LDOS of a corresponding microscopic model Hamiltonian, 
given by [lOt ] : 

H = 2 Yl %4r c i* + 2 ^2( A ij c n c n + cx -) - v Yl c ^° Ci ° ( 10 ) 



Here Yluj) denotes a sum over distinct pairs of nearest neighbor atomic sites on a square 
lattice with N sites, c] CT creates an electron with spin a (| or J,) at site j, /j, is the chemical 
potential, and Ay- denotes the strength of the pairing interaction between electrons at sites 
i and j. Finally, we write £jj as 

% = -^ op e * A 5, (11) 



with 



^=Y c f J Mrl-df. (12) 



Here the integral runs along the line from atomic site z, located at i, to the atomic site j, 
located at j, and t hop > is the hopping integral for nearest neighbor sites on the lattice. 
Note that the prefactor in A'^ involves the factor of hc/e due to a single electronic charge, 
and is thus twice as large as that in Ay, which involves the charge of a Cooper pair. 
Following Ref. [lo| we take A^ to be given by 

a,, = \ \M+iM e ^ 7 (is) 



where 



and 



9i + 0j)/2, if bond is in x-direction, 
% - { " (14) 

?j + 0j)/2 + 7T, if bond is in y-direction, 



Ai- A i ( ^ ; (15) 



is the value of the complex superconducting order parameter at site j. We will refer to 
the lattice over which the sums in fflOl) are carried out as the atomic lattice (in order to 
distinguish it from the XY lattice.) The first term in Eq. ( TlOl) thus corresponds to the 
kinetic energy, the second term is a BCS type of pairing interaction with d-wave symmetry, 
and the third is the energy associated with the chemical potential. 

The model we present is similar to the one presented in Ref. lOj, the main differences 
being the inclusion of a vector potential in the GL free energy functional, and the spatial 
distribution of the inhomogeneities. Because the vector potential introduces vortices in the 
system, our results differ substantially from our previous work. 
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III. RESULTS 



We first present results for homogeneous systems in the presence of an applied transverse 
magnetic field equal to two flux quanta through the lattice at low T. This is the lowest 
magnetic field consistent with the periodic boundary conditions. In this case, t c0i and 
Aj(0) are independent of i. Fig. [1] shows our calculated maps of A for a homogeneous 
system. In part (a), the lengths and directions of the arrows represent the magnitude and 
phase of A in each XY cell. Even though t c oi is homogeneous, the magnetic field renders 
A inhomogeneous, especially near the vortex cores, where A has a large phase gradient 
and a smaller magnitude. This behavior is familiar from Ginzburg-Landau treatments of 
homogeneous Type II superconductors in a magnetic field. Part (b) of Fig. [T] shows a map 
of |A|, with dark (light) regions representing small (large) value of |A|. The vortex cores 
are the darkest regions. 

Fig. [2] show the LDOS averaged over regions near and far from the vortex cores of the 
system described in Fig. [TJ Since |A| is small near the core, we have defined an XY cell to 
be "near the core" if |A| < |A avg |/2 in that cell, where |A avg | is the value of |A| averaged 
over all the XY cells of the lattice. All other cells are considered to be far from the core. 
Thus, in order to compute the averaged quantities shown in Fig. [2j we first calculated the 
LDOS on every atomic site by exact numerical diagonalization of the Hamiltonian ffTUl) . 
and then averaged the LDOS over the the set of atomic sites near and far from the vortex 
cores, as defined above. We have chosen to enclose nine atomic sites inside each XY cell for 
this system, because the coherence length in cuprates is ~ 15 A, approximately three times 
larger than the atomic lattice constant, ~ 5 A. 

Fig. [2] shows that, in a homogeneous system, the LDOS far from the core is strongly 
suppressed near u = 0, and exhibits sharp coherence peaks, reminiscent of a d-wave super- 
conductor in the absence of a magnetic field. Near the vortex cores, on the other hand, 
the gap is filled, and the LDOS is large near u = 0, resembling the spectrum of a gapless 
tight-binding model in two dimensions and zero magnetic field [10, |34j- The spectrum near 
the core has considerable numerical noise, because it represents an average over only a few 
atomic sites. In particular, the rapid oscillations are probably due to this numerical artifact. 

To estimate our magnetic fields, we note that, because of the numerical implementation 
of the periodic boundary conditions 10, [35J, the field has to be chosen so that the number 
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N v of magnetic flux quanta $o — he/ (2e) through the system is a multiple of two. Therefore, 
the magnetic field can be estimated using B = N v &q/S where S is the area of the system, 
and $o = he/ (2e) ~ 2 x 10~ 15 T-m 2 . In the present article we use a numerical sample of area 
area S = (48ao) 2 , where a ~ 5 A is the atomic lattice constant. Therefore B w N v x 3 T, 
and for the systems containing two vortices, B ~ 6T. The magnitude of this magnetic field 
is comparable used in STM experiments Q- 

We now briefly discuss the temperature evolution of (|A| 2 ), defined as an average of 
(|Aj| 2 ) over all XY lattice cells i, for systems with homogeneous t c o in a transverse magnetic 
field. Here, (...) denotes a thermal average, obtained using Monte Carlo simulations. Fig. [3] 
shows curves of (|A| 2 )(t) versus reduced temperature t, for systems of the same size subject 
to different magnetic fields, and therefore containing different numbers of vortices N v . In 
the zero magnetic field case, N v = 0, (|A| 2 )(t) has a minimum as a function of t, which 
occurs near the phase ordering temperature. Fig. [3] shows that, as the field is increased, the 
phase ordering temperature is reduced, and around N v = 36, it seems to drop to zero. The 
fact that of (|A| 2 )(t) increases with t for large t is an artifact of the GL functional, as has 
been discussed in detail in [10] for the case B = 0. 

We now proceed to show systems with inhomogeneities. Fig. H] shows gap maps of systems 
in which t c0i , and therefore JxY,ij, are are i dependent. Specifically, two square regions, of 
size 4x4 XY cells each, have t c0i = ft c0 , with / = 3; those are called (3 cells, as described 
in the previous section. The remainder of the XY cells have £ c oi = t c o, and are called a 
cells. Clearly, the vortices are pinned in the (3 regions, where |A| is large. The fact that |A| 
is large in f3 regions is of course due to the fact that at low temperatures, |A| is roughly 
proportional to t C Q. On the other hand, in order to understand why the vortices are pinned 
in the large-|A| regions, we note that JxY,ij at low temperatures can be estimated with the 
use of eq. ([9]). For the particular parameters used in this calculation, namely, x = 3 and 
/ = 3, JxY,ij is about 27 times smaller within (3 regions than in the a regions. Because 
of this ratio, phase gradients near the vortex cores cost much less energy in the f3 regions 
than in the a regions, even though |A| is larger in the f3 regions. Thus, it is energetically 
favorable for the vortices to be pinned in the j3 regions. 

Fig. [5] shows the LDOS averaged over a and (3 regions of the system shown in Fig. HI 
This Figure shows that, far from the vortex cores (a regions), the LDOS is very similar 
to that of regions far from the core in homogeneous system at low magnetic field [cf. Fig. 
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[2] - both spectra have a small gap and sharp coherence peaks. However, in contrast to 
the homogeneous case, the LDOS of the inhomogeneous system exhibits a large gap and 
broadened peaks near the vortex core. This LDOS spectrum is similar to that calculated 
for the large gap regions of systems with quenched disorder and no magnetic field in our 
previous work 10j. In the present case, the gap in the LDOS is large in the (3 regions, of 
course, because of the large value of t c oi and therefore of |A|, near the vortex cores. The 
peaks in the LDOS in the (3 region above the gap are low and broad, on the other hand, 
because of the large phase gradient of the superconducting order parameter in these regions. 
It is as if the system has lost phase coherence in those regions due to the presence of the 
vortex core. We can more clearly see this point by comparing the LDOS of a (3 region 
containing a pinned vortex core to that of a (3 region in a system in a zero magnetic field, 
in which, therefore, there is no vortex core to be pinned. We now describe this system. 

Fig. [6] shows a system at low temperature with two inhomogeneities, similar to that shown 
in Fig. HI but with zero rather than a finite magnetic field. Clearly, in the ground state, as 
expected, the phase of A is almost uniform. Furthermore, and also as expected, | A| is larger 
in the (3 regions, which have a larger value of t c o«. However, when we turn to the LDOS 
(Fig. [7]), the LDOS of the (3 regions is characterized by much sharper coherence peaks than 
that at finite magnetic field described in the previous paragraph, presumably because of the 
absence of a large phase gradient. 

We have also tested the sensitivity of our results to lattice size and to the number of 
atomic sites per XY cell. To do this, we performed calculations similar to the ones described 
above, but with 24x24 instead of 16 x 16 XY lattices, each XY cell containing four instead of 
nine atomic lattice sites. Fig.[8]shows the results for a 24x24 XY lattice of an inhomogeneous 
system in the presence of a magnetic field. As in the 16x16 case, the vortices are pinned 
in the regions with large A and small JxY,ij- Fig. [9] shows that the LDOS for this system 
both near and far from the core is nearly the same as that of the 16x16 system shown in 
Fig. [5j Likewise, a system with 24x24 XY cells and two atomic sites per XY cell, but with 
no magnetic field, is shown in Fig. [101 the corresponding LDOS averaged over regions near 
and far from the core is shown in Fig. [TTJ The LDOS shown in this Figure has the same 
features as that of the 16 x 16 system shown in Fig. [7J We thus conclude that our results are 
not very sensitive to the lattice sizes used, nor to the number of atomic sites per XY cell. 

In the calculations described above, have assumed that A 2 oc T c q. If y > 3, then the 
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vortices are consistently pinned in the large gap region. By contrast, if y = 1, we find the 
vortices may or may not be pinned in the large gap region, depending on where these regions 
are located within the computational lattice. The vortices obviously go to the large gap 
region because of the correlation between A 2 and T c0 . As mentioned earlier, this correlation 
originates in the fact large A 2 implies a small energy cost to introduce a gradient in the 
phase of the order parameter. Such a gradient must exist near a vortex core; therefore, the 
vortex prefers to be in a region where this gradient is energetically inexpensive. Melikyan 
and Tesanovic jsoj discuss various other possible causes of larger pinning in the large gap 
region; these include quantum fluctuations, distortion of the atomic lattice by the vortices, 
and the pinning of vortices by impurity atoms. Our model could be viewed as a special 
kind of such impurity pinning, in which the "impurities" are superconducting regions with 
a large gap and large penetration depth. We have also looked at how the size of the pinning 
regions affects the pinning. The pinning is generally more effective for large pinning regions, 
probably because the vortex energy is reduced by a larger amount if pinned in a region of 
large pinning area, all other parameters being the same. 

Finally, we note the similarity between our calculated LDOS versus energy curves for 
regions near the vortex cores, and the corresponding curves for large gap regions of inhomo- 
geneous systems with quenched disorder in a zero magnetic field, calculated in our previous 
work [lo|. Those similarities have been implied in several experimental papers. For example, 
Lang et al. J] discussed the similarity between the spectra of large-gap regions of inhomo- 
geneous systems at low temperatures in a zero field and those observed in the pseudogap 
regime of some cuprate superconductors. In the pseudogap region, the phase configuration 
is disordered and therefore, possibly like the j3 re gion at zero field, both might be consid- 
ered as "normal" regions. Likewise, Fischer et al. [36( have also noted the similarity of the 
low-T spectra near the vortex cores to spectra in the pseudogap regime. Once again, the 
similarity may arise because both the interior of the vortex cores and the pseudogap region 
may be considered as "normal" . Thus, both of these reports implicitly suggest a similarity 
between the low-T spectra in the large gap regions of a disordered system at B = and the 
corresponding spectra near vortex cores at finite B. 
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IV. DISCUSSION 



We have presented a model for the pinning of vortices by inhomogeneities in a two- 
dimensional type II superconductor at low temperatures. The model is based on a GL 
free energy functional, and it is inspired by our previous study jlo| of the LDOS in an 
inhomogeneous superconductor in zero magnetic field. In our model, we have proposed a 
GL free energy functional in which regions with a large value of the superconducting order 
parameter have a large penetration depth, as suggested by zero-field STM experiments on 
cuprates. Using an annealing process and Monte Carlo simulations, we have found that the 
pinning of vortices by those large-gap regions emerges naturally from the functional form 
of our GL free energy. The vortices are attracted to the large-gap regions, in our model, 
because the large penetration depth leads to a small coupling between cells in those regions. 
Therefore, the phase of the superconducting order parameter can more easily bend in these 
regions, and this in turn allows the large-gap regions to accomodate a vortex more easily 
than regions with a small gap. By contrast, in the absence of quenched inhomogeneities, 
minimization of the GL free energy functional yields a spatial configuration in which the 
superconducting order parameter is suppressed near vortex cores. 

It is worth commenting further on the qualitative physics underlying the pinning of the 
vortices in the region of enhanced local T c . Basically, this pinning behavior occurs because, in 
our model, a locally enhanced Tc corresponds to a locally suppressed superfluid density. This 
superfluid density is proportional to the local 1/A 2 , and is related to the Ginzburg-Landau 
coefficients, as we have described earlier. The vortices can be more easily accommodated 
in these large-T c regions because the large phase gradients which characterize the vortices 
cost less energy in such regions. Although we describe the pinning in terms of the variation 
of 1/A 2 , the pinning is not due to any kind of "magnetic" forces - the quantity 1/A 2 is 
just a way of describing the local superfluid density. Thus, just as in conventional pinning, 
the vortices are attracted to regions of lower free energy. The major difference is that the 
superfluid density (and 1/A 2 ) is proportional to T c in the present model, rather than being 
independent of, or inversely correlated with T c as in more conventional pinning. This leads 
to pinning in a large-T c region rather than a small-T c region. We also emphasize that, 
although our pinning results are obtained using rather elaborate numerical calculations, the 
underlying physics is straightforward: the pinning behavior is consistent with qualitative 
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expectations, given the model. 

It may seem strange to consider the spatial variation in 1/A 2 as occurring over a scale 
of £o ~ 15 A, when A itself is of order 10 3 A or more. However, it should be remembered 
that 1/A 2 is related to the coefficients of the Ginzburg-Landau free energy [see eq. (1)], 
and these coefficients are, in fact, expected to vary over a scale of £o- Thus, the model 
is, indeed, reasonable and consistent with the expected physics. It is the local superfluid 
density (proportional to 1/A 2 ) which varies over a scale of £o- 

It is also worth commenting on what kind of real systems could be described by our 
model. Even if one neglects the weak Josephson coupling between the layers of a high- 
Tc material, there will still be magnetic interactions between pancake vortices in adjacent 
layers, so that the system will not be truly 2D. Likewise, if we wish to use this model to 
describe a very thin 2D high-T c film, there are stray fields in vacuum extending into the third 
dimension, which will contribute to the total energy. In both cases, our model is definitely an 
oversimplification. Nonetheless, our model does seem to describe some observed features in 
real high-T c materials, suggesting that it captures some significant physics in these systems. 
Thus, we consider our model as a possible starting point for a fully realistic treatment of 
either 3D high-T c materials or very thin 2D high-T c layers. 

We have connected our work on the pinning of vortices to continuing efforts by several 
groups to describe the density of states near vortices in cuprates superconductors. To do this, 
we have studied a model BCS Hamiltonian with d-wave symmetry, in which the pairing field 
is obtained from simulations of the GL free energy functional. We use exact diagonalization 
to compute the local density of states on each atomic site of the lattice described by this 
Hamiltonian. For homogeneous systems, we found that the LDOS near the vortex cores 
resembles that of a gapless tight-binding model in two dimensions, with a Van Hove peak at 
zero energy. However, when we introduce the inhomogeneities with a large pairing field and 
large penetration depth that pin the vortices, the LDOS near the vortex cores is markedly 
different from that of the homogeneous systems. Namely, this LDOS exhibits a large gap 
as well as low and broadened peaks at energies greater than that of the superconducting 
gap fa.r from the vortex core. Also, our calculated LDOS near the vortex cores in this 
inhomogeneous case does not exhibit the spurious zero-energy peak which is present in 
several other theoretical studies but is absent from experiments. All of these features in our 
calculations are consistent with results observed in STM experiments on YBCO JjJ, ll2j and 
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B Q, H, 



Bi2212 [13|,l]J,ll5!,ll6|. 



30t ] , using a somewhat different 



Our results are also consistent with those obtained of Ref. 
model. Those authors obtain better agreement between their calculated LDOS spectra and 

n 

experiment [14[] , if they assume an enhanced rather than uniform pairing field near the vortex 
cores. In particular, introducing this large pairing field near the vortex core in their model 
suppresses the unphysical zero-energy peak. Our approach, in which a large gap is correlated 
with a large penetration depth, may heop justify the occurrence of this larger pairing field 
in the vortex cores. 

Finally, we briefly comment on the possible physical origin of correlation between large 
gap and large penetration depth. The origin of the spatial fluctuations of T c0 observed in 
the cuprates may be the local fluctuations in concentration, which are highly likely since the 
cuprates are mostly disordered alloys and the local concentration would involve an average 
a coherence length, which is only about 15 A . But T c0 is proportional to the gap (i. 
e., presumably, the pseudogap), which increases with decreasing concentration of charge 
carriers, whereas 1/A 2 (0) is proportional to the superfluid density, which should decrease 
with decreasing charge carrier concentration. Therefore, we expect T c q to be positively 
correlated with A 2 (0), as seen experimentally and as used in the present model. 
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FIG. 1: Map of the pairing field A over a homogeneous, two-dimensional superconductor of 
area 16£o x 16£o m a transverse, uniform magnetic field B ~ 6T. The superconductor has been 
discretized into XY cells, each of which has an area £q and encloses nine atomic sites. In part (a) 
the length and direction of each arrow represents the complex value of A within an XY cell. Part 
(b) shows a map of the magnitude |A| of the pairing field. Dark (light) regions represent a small 
(large) value of |A|. Vortex cores can be easily identified as the darkest regions. 
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FIG. 2: Local density of states averaged over two different groups of XY cells in the homogeneous 
superconductor described in Fig. [1} far from the vortex core (full curve) and near the core (dotted 
curve). The oscillations in the dotted curve are due to numerical fluctuations arising from the 
small number of sites in the core regions 
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FIG. 3: Thermal average (|A| 2 ) of the squared magnitude of the superconducting order parameter 
versus the reduced temperature t for a system with a homogeneous t C Q placed in various fields. The 
magnitude of the magnetic field is B ~ N v x 3T as described in the text. We can observe that the 
minimum of (|A| 2 )(t) versus t (which is known to occur near the phase ordering temperature in 
zero magnetic field systems [10I.I37I]) is shifted toward smaller values of t with increasing magnetic 



field. 
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FIG. 4: Map of the pairing field A in an inhomogeneous, two-dimensional superconductor of area 
16£o x 16£o m a transverse, uniform magnetic field B ~ 6 T. The superconductor has two (3 regions 
where t C Qi is large; these correspond to the light regions in (b). The system has been discretized 
into XY cells, each of which has an area £q and encloses nine atomic sites. In part (a) each arrow 
represents the complex value of A within an XY cell. Part (b) shows a map of the magnitude |A| 
of the pairing field. Dark (light) regions represent a small (large) value of |A|. The locations of 
the two vortex cores can be identified as the regions with a large phase gradient in part (a). They 
are pinned to regions with a large |A| because of the low value of the coupling between XY cells 
in those regions. The present results are obtained from Eqs. ([8]) and Q, using / = 3 and x = 3. 
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FIG. 5: Local density of states averaged over two different groups of XY cells in the inhomogeneous 
superconductor shown in Fig. [H far from the core (a) and near the core (/?). In agreement with 
experimental results, this Figure shows that (i) far from vortex cores the LDOS shows sharp 
coherence peaks; (ii) near the vortex cores, the LDOS does not display the unphysical zero-energy 
peak obtained in other models. Instead (iii), it has a large gap, as well as low and broad peaks. 
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FIG. 6: Same as Fig. U but for a system at B = instead of B ^ 0. At low temperatures, the 
phase of A is almost uniform, since, in the absence of a magnetic field, phase gradients cost energy. 
As expected, A is larger in the (3 regions, which have a large value of £ c oi- 
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FIG. 7: Local density of states for the system shown in Fig. [6l The LDOS in the (3 regions of this 
system has much sharper peaks than the LDOS of the j3 regions for the corresponding system in a 
finite magnetic field shown in Fig. [5j This result shows that having a large gap is not a sufficient 
condition to observe broadened peaks near vortex cores; a large phase gradient is also required. 
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FIG. 8: Similar to Fig. HJ but the system has 24x24 instead of 16x16 XY cells, each with four 
instead of nine atomic lattice sites. 
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FIG. 9: Local density of states averaged over regions far from the vortex cores, for the system 
shown Fig. [HJ Results are very similar to the corresponding system with 16x16 cells, which has 
nine instead of four atoms per XY cell, and is shown in Fig. [5j Thus, our results are not strongly 
dependent on the size of the XY cell or of the XY lattice. 
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FIG. 10: Same as Fig. [6] but for a system consisting of 24x24 instead of 16x16 XY cells, each 
having four instead of nine atomic sites. 
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FIG. 11: Same as Fig. [7] but for a system of 24x24 instead of 16x16 XY cells, each having 
has four rather than the nine atomic sites. 
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